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We propose a theory for a type of apertureless scanning near field microscopy that is intended to 
allow the measurement of magnetism on a nanometer length scale. A scanning probe, for example 
a scanning tunneling microscope (STM) tip, is used to scan a magnetic substrate while a laser is 
focused on it. The electric field between the tip and substrate is enhanced in such a way that the 
circular polarization due to the Kerr effect, which is normally of order 0.1% is increased by up to 
two orders of magnitude for the case of a Ag or W tip and an Fe sample. Apart from this there is a 
large background of circular polarization which is non-magnetic in origin. This circular polarization 
is produced by light scattered from the STM tip and substrate. A detailed retarded calculation for 
this light-in-light-out experiment is presented. 

61.16.Ch, 78.20.Ls 



I. BACKGROUND 

A microscope that can measure magnetic structure 
with high resolution would be very desirable for practical 
reasons such as investigating read heads and media for 
magnetic recording and because magnetic nanostructures 
are currently of great interest. We propose a scanning 
magnetic micrcficope with a resolution of approximately 
1 nm to 10 mm depending on the tip shape. This device 
is a form of apertureless near-field scanning magneto- 
optical microscope in which light from a laser is focused 
at the tip of a scanning tunneling microscope (STM) or 
atomic force microscope (AFM) or at an isolated sphere 
on the sample surface and the circular polarization of 
the light scattered from the tip and sample is measured. 
This is analogous to attempts to use the circular polar- 
ization of the light emitted by a scanning turuoeling mi- 
croscope to construct a magnetic microscope.ETH It was 
found that the light-emission scheme is not suitable for 
use as a magnetic microscope due to the small degree of 
circular polarization produced by a magnetic sample as. 
well as the relatively low intensity of the emitted light .□ 
The light intensity is low despite the fact that the electric 
field at the tip of an STM can be strongly enhanced by 
interface plasmons localized between the sharp tip apex 
and the substrate. Without the enhancement the light 
intensity would be smaller by several orders of magni- 
tude. Other techniques to develop magnetic microscopes 
based on apertureless near-field saaiuiing have resolution 
on the order of at least 100 nm.Bi3 Attempts to mea- 
sure magnetic susface structure by means of the STM 
tunneling currcntB have also been made and recently the 



magnetic structure of an antifesromagnetic overlayer was 
mapped using this technique.ta 

In the scheme we propose, the scattered laser light in- 
tensity is large compared to that produced in STM light 
emission experiments making detection far easier. The 
circular polarization of the scattered light has both a 
component that is independent of the state of magne- 
tization of the sample and a component that depends on 
the magnetization. The latter effect is rather small both 
when light is scattered from a magnetic surface as in the 
surface magneto-optic Kerr effect, and when light is emit- 
ted from an STM. However, when a tip is present and a 
laser provides incident light that is linearly polarized par- 
allel to the surface (s polarized) and propagating in the 
direction of magnetization, we find up to two orders of 
magnitude increase in the degree of circular polarization 
compared with the basic Kerr effect which is generally 
of order 0.1 %. This enhancement occurs because the s 
polarized light that scatters from the tip and substrate 
develops a component in the direction perpendicular to 
the surface that is out of phase with the s polarized light. 
This is due to (1) the ordinary Kerr effect in which some 
of the incoming s polarized light is converted to p polar- 
ization and (2) the Kerr effect that acts on the light that 
is scattered between the tip and the substrate in the near 
field. The components of these fields that are perpendic- 
ular to the surface and lie between the tip and sample 
are enhanced by one to two orders of magnitude with a 
corresponding increase in the circular polarization. 

The part of the circular polarization that does not de- 
pend on the magnetism of the surface is due to light that 
is reflected from both the tip and sample surface. The 
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surface reflection coefficients for s polarized and p polar- 
ized light result in phase differences that produce circu- 
lar polarized light. This circular polarization is large and 
has a characteristic dependence on the observation angle 
which is different from the one displayed by the circular 
polarization that is due to magnetic effects. This differ- 
ence in angular dependence may help in separating the 
two contributions. 

The resolution of the microscope is roughly set by the 
radius of the tip (probe). A smaller probe gives a bet- 
ter resolution, however, one must also take into account 
that a very small probe gives a very small scattering cross 
section. At the same time, we find that a small probe 
is essential in obtaining a high degree of magnetism- 
dependent polarization of the scattered light. Our results 
indicate that a probe size of about 20 nm is optimal in 
striking the balance between these conflicting require- 
ments. 

The paper is organized in the following way. In section 
U we present a general theory with which we can solve 
the multiple-scattering problem defined by the model tip- 
sample geometry and calculate the scattered radiation 
and its polarization state. This theory uses a spherical 
model tip. The electromagnetic field is calculated tak- 
ing retardation effects fully into account, while the mag- 
netooptic corrections to the reflection coefficients of the 
sample are treated to flrst order in the off-diagonal ele- 
ments of the dielectric tensor. In section III, we discuss 



the limit of a point dipole tip, in order to demonstrate the 
central physics more clearly, and also to make a compari- 
son with the case ofJight emission from an STM probing 
a magnetic surfacccl In Sec. |^ we present and discuss 
the numerical results obtained using the general theory 
of Sec. 0. 



II. THEORY 
A. Basic considerations 

In this section we present a theoretical formalism that 
allows us to calculate the polarization of light scattered 
from a tip-sample geometry in which the model tip is a 
sphere of finite size (the radius of the sphere corresponds 
roughly to the radius of curvature of the real tip). The 
calculation takes into account higher multipoles on the 
sphere as well as retardation effects. 

We consider a situation illustrated in Fig. [| where the 
magnetic sample fills the half-space 2: > 0, the model tip 
is a sphere with radius R that has its center on the z axis 
at 2: = — (rf+i?), thus the tip-sample separation (smallest 
distance) is d. A plane wave incident on this system with 
wave vector q can be written as (all fields etc. vary with 
time as e~*"*, but we omit this factor here and in the 
following) 



E(r) ^ |£;(")[z X q||] -£;(p)[qx (z x q||)]|: 



,iqr 



(2.1) 



For any given wave vector q the polarization vectors for 
s and p polarization are, respectively. 



and 



s = z X qii 



P = -q X (z X q||). 



where qy is the projection of q in the surface plane. Both 
s and p are orthogonal to the unit vector q. We will cal- 
culate the radiation sent out in an arbitrary direction q' 
after refiection off the tip-sample system. The radiated 
electric field propagating along q' can be written as 



Erad — 



{<1[^xq[|]-i?(r^[q'x(zxqj|)]}. 



(2.2) 



The polarization state of this radiation is different from 
that of the incoming light. First, there is a rather large, 
in this context less interesting, contribution to the circu- 
lar polarization that comes about because the reflection 
probability from the surface is different for s polarized 
and p polarized light. In addition, the scattered light 
acquires a certain degree of circular polarization due to 
the Kerr effect; of course, for our purposes this is the in- 
teresting contribution. An incident s polarized wave will 
electrically polarize the tip in a direction parallel to the 
sample surface, but reflections off the magnetic sample 
surface also gives rise to a fleld component perpendicu- 
lar to the surface. There is a strong enhancement of the 
component of the fleld perpendicular to the surface due 
to the cavity formed between tip and sample, and there- 
fore this second contribution to the circular polarization 
should also be detectable. 

The degree of circular polarization of the outgoing ra- 
diation is given in terms of the Stokes parameters as 



^3 

PCP - TT- 
Do 



(2.3) 



Following Jackson's deflnitionsjl^ the Stokes parameter 
S'o is proportional to the total radiated differential power, 
whereas 6*3 is proportional to the difference in intensity 
between left and right hand circularly polarized light. 
With the geometry illustrated in Fig. || one finds that in 
the radiation zone, far away from the tip and the sample 
surface, s = (j)' and p = 9', where cj)' and 6' denote the 
angular direction of the scattered light. 

Left hand circularly polarized light (positive helicity) 
is associated with the unit polarization vector e+ — 
{p+is)/y/2, while for right hand circularly polarized light 
(negative helicity) e_ = (p — is) / \/2. For the radiation 
propagating along the q'-direction 
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Hence the degree of polarization can be expressed as 



pcp = -2Im 



rad rad 



rad I J 



(2.6) 



B. Sample dielectric tensor 

We have expressed the measured magnetic circular 
dichroism in terms of the field amplitudes at the detec- 
tor. These in turn have to be calculated with the tip- 
sample aeametry and tip and sample dielectric functions 
as input.ll3 Because the sample is magnetic, its dielectric 
function is a matrix that takes the form 



es(t^) 
— ei cos ip 
ei sin 7 sin ip 



ei cos Lp 
es(w) 
-ei cos 7 sin 93 



— ei sin7sm(/3 
ei cos 7 sin 
es(^) 



(2.7) 



where es(^) is the substrate dielectric function. (Some- 
times ei appearing in the off-diagonal elements is writ- 
ten ei{oj) = iQ{uj)es{L^j), Q being the so called magneto- 
optical constant.) The angles (p and 7 relative to the x 
axis specify the direction of the magnetization. We will 
work in a configuration where 1^9 = 7 = 7r/2 which means 
that the magnetization is directed along the y-axis 

M = MM = My, 

(see Fig. |l|). Tiien the dielectric tensor for the sample 
rnM 



takes the forr 



esiuj) -61(0;) 

es{co) 
ei(w) es{Lo) 



(2., 



In Table || we list values, taken from Ref. |T^, for ei{uj) 
used in the present calculations. 



C. Surface response 



When a plane wave of the form in Eq. (2.1) is reflected 
from the sample surface, q|| remains unchanged. The 
reflected wave can in general be written as 

E, = {i?(^)[z X q||] - i?(f)[q x (z x q,,)]} e'^^--, (2.9) 

where 

'ill, 



(qih-p) with ./k'^-qf 



and k = u/c, the subscript — expressing that the re- 
flected wave propagates in the negative z direction. 

The reflected amplitude varies linearly with the ampli- 
tudes of the incident wave. Thus, for s polarization, 



p'spE^"^- 



(2.10) 



The response function ps describes the "normal" reflec- 
tion that takes place whether or not the sample is mag- 
netic. The Fresnel formula for s polarization is 



P-PS 
P + Ps' 



(2.11) 



where 



PS 
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k = io/c, and ks = ky^es{uj) are the wave- vector magni- 
tudes in vacuum and the sample material, respectively. 
The second term in Eq. (2.10) applies to magnetic sur- 
faces for which the off-diagonal element ei(w) in the di- 
electric tensor is non-zero and polarization mixing as a 
result of "anomalous" reflection events becomes possible. 
The conversion factor can be explicitly written as 



ei(w) q\ipk 



sp 



PS ip + Ps){^s{^)p + Ps, 



■[M-q||] = pL[M-q||]. 

(2.12) 



This expression caft,be derived by the methods devel- 
oped by Zak et o/.li3 There is also a contribution from 
the off-diagonal elements in the dielectric tensor to the 
s-to-s reflection coefficient. However, this contribution is 
exceedingly small since it depends on the square of ei{uj) 
and we neglect it in the follow ing. 

In analogy with Eq. ( 2.10 ), the reflected p polarized 
light is associated with the electric field 



EP = p,E^P> + p'^,E^'^> + p;^E^P\ 



In this case the Fresnel formula is 
es{uj)p-ps 



Pp 



es{uj)p + ps' 



(2.13) 



(2.14) 



while the conversion factors originating from £1(0;), p'p^ 
and p'pp are given by 



P'ps = -P'sp = -Pl[M ■ q||] 
[pL was introduced in Eq. ( ^.12 )] and 
, 2ei{uj)pq 



(2.15) 



{€s{uj)p + ps) 



[z-(q|| xM)] = pT[z-(q|| xM)]. 

(2.16) 



Thus, unlike p'^^ and p'^^, p'^^ vanishes when qy and M are 
parallel. Figure ^illustrates the three possible anomalous 
scattering processes. 
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The expressions for all the surface response functions 
can be extended to the case of evanescent waves for which 
(7|| = |q[| I > fc, and p is imaginary. (When evaluating the 
square roots defining p and ps the branch cut lies below 
the positive real axis.) Later when we deal with the cou- 
pling between the sphere and the sample this becomes 
important since the coupling is mediated both by prop- 
agating and evanescent waves. 



D. Scattering from the sphere 

We now deal with the scattering off a sphere whose 
optical properties are described by its dielectric function 
erii^)- To this end we expand the electromagnetic field 
inside, and just outside the sphere in terms of electric (E) 
and magnetic (M) multipoles. We use the same notation 
as Jackson,EJ so that inside the sphere the electric field 
is written (the spherical coordinate system used here has 
the origin at the center of the sphere) 



Im 



-V X 



-■Im 



ji{kTr)lLi, 



(2.17) 



while just outside the sphere the solution is a linear com- 
bination of incident and outgoing waves which can be 
written 



E 



E 

Irn 



V X [[a\^n{kr) + b\^hi{kr)] X,„} . (2.18) 
The corresponding magnetic field outside the sphere is 



B 



-V X 



{ ^jK^-^) + b^w.?hi{kr)\ X,„,} . (2.19) 



Here the vector spherical harmonics are defined as 



x,™(0,0) = [Ly,™(0,0)]/yi(7TT), 



(2.20) 



L = — ir X V is the angular momentum operator, and Yim 
are the usual spherical harmonics. The vector spherical 
harmonics form an orthonormal set of functions on the 
surface of a unit sphere. The functions j; denote spheri- 
cal Bessel functions, while hi = /i^' are spherical Hankel 
functions describing outgoing waves. As before k = lo/c, 
kx — v^erfc, and r denotes the distance from the cen- 
ter of the sphere. The coefficients cim, aim, and bim for 
the different multipoles are as yet unknown; they will be 
calculated in the following. 

Since the angular momentum operator L does not have 
any vector component in the radial direction, the same 



holds true for X;™. Consequently, for an electric mul- 
tipole field, the magnetic field does not have a radial 
component. On the other hand the electric field, as one 
should expect by visualizing the field from an electric 
dipole, has a radial component. For a magnetic multi- 
pole field, the roles are interchanged; the magnetic field 
but not the electric field has a radial component. 

In the present context, we are primarily interested in 
knowing how the sphere acts as a scatterer. It therefore 
suffices to know what outgoing waves are obtained for 
given incident waves. The ratio between the b and a co- 
efficient for each multipole provides this information. By 
demanding that the tangential E and H fields and the 
normal B and D fields are continuous we arrive at the 
following expressions for the sphere "response" functions 



er kR][{kR)+3i{kR) {er ~ I ~ Ji) 



er kRh[{kR) + hi{kR) {er - 1 



and 



AM) ^ Ol„r _ 



kRj[{kR)-ji{kR)Ji 
kRh[{kR) - hi{kR)Ji' 



-JiY 
(2.21) 

(2.22) 



where Ji is shorthand for 

Ji = kTRj'i{kTR)/ji{kTR). 

Thanks to the symmetry of the sphere these response 
functions are independent of m. 

One additional ingredient is needed before we can cal- 
culate the field around the sphere when it is interacting 
with the sample. We have to be able to calculate "over- 
lap integrals" between plane waves and multipole fields 
and vice versa. (Note that for lack of a better terminol- 
ogy, we will use "plane wave" as a common name both 
for propagating plane waves and waves that propagate 
in the directions parallel to the sample surface but are 
evanescent in t he t hird, z direction). A plane wave, like 
the one in Eq. (2.1), impinging on the sphere can be ex- 
panded in terms of r egular mul tipole contributions [the 
ji terms in Eqs. (2.18) and (2.19)]. The resulting a coeffi- 
cients depend linearly on the amplitudes of the incoming 
wave 



a 



(M) 



and 



is) 



(2.23) 



Explicit expressions for the proportionality factors / are 
derived in the Appendix. 

We also have to go in the opposite direction, from a 
certain multipole to a plane wave. In general the field 
radiated from the sphere can be written as 



E(r) 



gf„-,(Q)(zxQ||)-C(Q)(Qx(^xQ||)) . (2.24) 
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Here Q = Q|| ± z ^k"^ — (5| (the sign is chosen accord- 
ing to the direction of propagation or exponential decay 
relative to the z axis) and a denotes a polarization [(E) 
or (M)]. The factors gi^iQ) and gf^(Q) are the contri- 
butions of a given multipole on the sphere, (/, m, cr), to 
the amplitude of a plane wave with wave vector Q. The 
explicit calculation of these factors is deferred to the Ap- 
pendix. 



E. Solution of the multiple-scattering problem 

We are now in a position to solve the multiple- 
scattering prohlem by using the same method as in a 
previous paper .EJ We solve for the fields at the surface of 
the sphere and then calculate the electromagnetic fields 
elsewhere, in particular the radiated fields found far from 
the tip and sample. 

In order to carry out the calculations, it is convenient 
to collect the a and b coefficients of Eqs. ( p. 18 ) and ( 2.19| ) 
into vectors a and b (in "multipole space" ) with the struc- 
ture 



■10 I ■ 



JM). 
^1-1 1 ■ 



(2.25) 



and similarly for b. At the same time a diagonal tensor 
s*, ca n be forme d fro m the sphere response functions in 
Eqs. ( 2.21 ) and ( 2.22 ), so that one can write 



b^ 



(2.26) 



The magnetooptic surface response functions p'p^, p'^^, 
and are small, and one can therefore treat the mag- 
netooptic effects by means of a series expansion in ei . We 
begin by calculating the fields to zeroth order in the mag- 
netooptic response. The field impinging on the sphere is a 
sum of three contributions (Figure 3) which in our vector 
notation can be written 



+ Na. (2.27) 

Here a'^" represents the amplitude of the field from the 
original incident wave; a"'"^ gives the amplitude from the 
incident wave reflected once from the sample; and finally, 
Nd represent waves sent out from the sphere itself that 
return to it after being reflected from the sample. Note 
that since the last term contains the exact field (a) it 
accounts for all multiple scattering events in which light 
is scattered between the sphere and sample an arbitrary 
number of times. 

With an incident wave given by Eq . (|2.lD , the elements 
entering a'^''' are found by using Eq. (|2.23D together with 



Eqs. (A6) and (A7). In an analogous way, the elements 
of a^«^ 



(B):rcf 



CM-)PpE^'^ + fin:{^-)PsE^'\ and 



0^'"" - /,t?(q-)Ppi?(^) + f,Z%ci-)psE('\ (2.28) 



Here PpE'-P'^ and psE^'^ are the amplitudes of the re- 
flected wave. The / functions now have q_ as an argu- 
ment since the direction of propagation relative to the 
z axis changes upon reflection. In the last term of Eq. 
( p. 27 ), is a tensor in multipole space. Its elements de- 
scribe how a wave with angular momentum I' and m' and 
polarization a' hits the sphere, is reflected, and returns 
to the sphere again after reflection from the plane, now 
with angular momentum I and m and polarization a. As 
a consequence, a tensor element is written 



Ira.l''. 



(2.29) 

In this equation a" stands for either s or p polarization, 
and 



Q± = Qii =t zi //c 



IQii 



Equation (2.29) can be derived form ally along the lines 
of the calculation presented in Ref. |l^, but it can also 
be understood in a more intuitive way as follows (see 
Fig. H for a schematic illustration) . The first factor s^, , 
connects the a coefficient of the initially incident mul- 
tipole to the corresponding b coefficient; then the func- 
tion 5(Q+) describes the overlap between the multipole 
I' , to', ct' and a plane wave, p the reflection of that plane 
wave off the sample, and /(Q_) gives the overlap be- 
tween the reflected wave and the multipole I, to, ct. Be- 
cause the propagation to the sample surface and back 
involves all wave vectors, a two-dimensional wave-vector 



integration is required. To evaluate Eq. (2.29) it is best 
to use cylindrical coordinates in which the angular inte- 
gration is trivial. Thanks to the cylindrical symmetry of 
the model geometry only tensor elements with m — m' 
are non-zero. The remaining |Qj|| integration runs from 
to oo so that the sphere-sample coupling is mediated 
by both propagating and evanescent waves. 



Equation (2.27) is solved by a matrix inversion (of 
course then only a finite number of multipoles can be 
retained) 



^1 -1 

a= 1 -N 



a 



-<rof\ 



(2.30) 



Knowing the coefficients aim we can calculate the elec- 
tromagnetic field everywhere to zeroth order in the mag- 
netooptic response. The various ingredients entering 
the calculation of the scattering processes are shown 
schematically in Fig. |[ 

We next proceed to calculate the corrections to the 
fields to first order in the magnetooptic response. To 
this end we introduce another set of coefficients aj|^'' and 



'■Im 



that vary linearly with ei, the off-diagonal elements 
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of the sample dielectric tensor. Together a'l^^ and a'l^^ 
form the vector a! in multipole space. One can solve for 
a' from the equation 



a = a + i\a , 



(2.31) 



in which all terms are first order in ei. The first term 
on the right hand side acts as a source term for the first- 
order-in-ei fields, while the last term accounts for multi- 
ple, normal (i.e. not involving the Kerr effect) reflections 
between the tip and sample; the fields appearing there, 
described by a', have already undergone an anomalous 
reflection. 

There are two contributions to d'°^^, one due to an 
anomalous reflection of the incident wave the very first 
time it hits the sample, and one due to waves incident 
from the sphere that are anomalously reflected, thus 



r'/iof 



Kd. 



(2.32) 



Both terms on the right hand side are first order in ei, 
the first one being the ordinary Kerr term and the seco nd 
the tip-induced Kerr term. In analogy with Eq. (2.28), 



/(£;),rcf_ rEp 



/,^p(q_)(p;pi?(rf-fp;,i?(^)) 



(2.33) 



and 



+flZ'i^^)p'.pE<^'^. (2.34) 

In these expressions we have suppressed the q|| depen- 
dence of Ppp, Ppg, and p'^p. This dependence becomes 

important when evaluating the tensor elements in K, 



^ II £aa 

(27r)2 



E 



(Q_)p;,,,,„(Q||)5,%„? (Q^ 



(2.35) 



Here a" and a'" stand for either s or p polarization (and 
to first order in ei). The v arious factors play 



the same roles here as in Eq. ( 2.29[ ) except, of course, 
for the fact that we now collect the effects of anoma- 
lous reflection events in the sample described by p'^„^n,. 
As in the previous case the angular integrations can be 
done analytically, and here, due to the angular depen- 
dence of p'^ ^1 , an element in K is nonzero if, and only if, 
m = m' ± 1. Once K has been calculated, we evaluate 



and then solve for a' 



d' = I - N 



= 1 - N 



?'rcf 



Kd 



(2.36) 



F. Radiated field 

Finally, we calculate the radiated fields, and the degree 
of polarization. To this end we add together all contri- 
butions to the radiated field that have scattered from the 
sphere at one time or another. 

To zeroth order in ei, the field is foimd by considering 
the outgoing waves coming from the sphere as well as 
the corresponding waves reflected (normal reflection) off 
the sample surface. The specularly reflected wave that 
has not interacted with the sphere is ignored since an 
experimental measurement that is sensitive to polariza- 
tion changes induced by the Kerr effect must avoid this 
otherwise dominating, direct contribution. 

The radiated field is calculated to first order in ei along 
the same lines using the aj^ terms instead of the aim 
terms in the multipole expansion as a source for the ra- 
diation. There is also an additional, in practice rather 
small, contribution to the first-order fields that results 
from anomalous reflection of zeroth-order waves coming 
from the sphere and scattering off the sample surface a 
last time before going out to infinity. 

The field radiated directly from the sphere into a direc- 
tion defined by the angles 9' and (j)' [{9' , (p') = ft'] can be 
calculated by means of a stationary phase approximation 
yielding 



g*'^'" ik cos 9' 



dir — 



2n 



Ima 



X ^'C(q-) + '^'.9^^(q-; 



(2.37) 



Here 



q'_ = /sq_ = fc sin 0' (x cos 0' + y sini/)') — k\ cos^'lz. 

(with an analogous definition, having a positive z com- 
ponent, for q'|_ used below), and 



s,(J\/) 

'■Im 



(M) , i{M) 



'■hn 



and 



'■Im 



''Im ' 



SO that the expression contains contributions to both ze- 
roth and first order in ei. The radiation resulting from 
normal or anomalous reflections from the sample sur- 
face can be calculated in a way analogous to Eq. (2.37). 
The only modifications are that one must (i) multiply by 
the appropriate reflection factor, and (ii) use q'^ instead 
of q'_ as an argument in g since the light in this case 
first propagates towards the sample after having left the 
sphere. 

These contributions can now be added together so that 
the polarization pcp , of the outgoing light that has been 
scattered ©ff the sphere at least once can be calculated. 
We obtainlia 



e: 



is) 

ad 



2tt 



E 



Imcr 

X [5r^(q'_)+Ps(fcsin0').9?^(q'+) 
-fsinfp;,(qj|)C 



iq'+) 



(2.38) 
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and 



e: 



(p) 

ad 



27r 



Irna 



x[C(q-) + Pp(fcsin^')C(q;) 

+ cos0V;p(ql|)C(q;) + sin0V;s(q'|)5f™(q;)] , (2.39) 

where the terms are organized according to where the 
last scattering event takes place. The first terms de- 
scribe waves that come directly from the sphere, the sec- 
ond terms give the contributions from waves that come 
from a final, normal scattering event from the sample, 
while the remaining terms originate from a final anoma- 
lous scattering event from the sample. The degree of 
circular polarizatio n is then found by inserting these ex- 
pressions into Eq. 



III. DIPOLE LIMIT 

For illustrative purposes it is useful to study the the- 
ory developed in the previous section in a limit where 
both the sphere radius and the sphere-sample distance 
are much smaller than the wavelength of light. When 
these conditions are fulfilled the sphere can be treated as 
a point dipole and re tarda tion e ffects are negligible. 

Returning to Eqs. ( 2.21 ) and ( 2.22| ) we find the follow- 
ing limiting behavior for the sphere response functions 
when both fci? < 1 and \kTR\ < 1, 



Ct — 1 



(2;-l)!!(2Z + l)!! ler + il + iy 



and 



i{kR) 



21+d, 



(2; + l)!!(2/ + 3)!!' 



(3.1) 



(3.2) 



Thus the behavior of a small sphere is dominated by its 
electric-dipole response s^^^ ^ . Therefore, in the rest 
of this section we will only discuss the electric dipole 
excitations of the sphere and the vectors a and h can 
effectively be reduced to just three components. 

We also need limiting expressions for the coupling ma- 
trix elements iVj^^p, TVj^^i, and iV^L^.i^i (the other 

remaining elements of TV vanish). An explicit evalu- 
ation of N^i^ii in which the large qy limits of pp — 
{es — l)/(es + 1), Ps — > 0, as well as / and g are used 
yields 



Nil ^ ^iMi 



EE 

1-1,1-1 



1 es-1 



i?3 



eT + 2 es + 1 8(rf + i?)3 



(3.3) 



This expression can be obtained from more intuitive rea- 
soning. If an electric field E — i?x polarizes the sphere, 
a dipole moment 



p = Ane.R^ 



(3.4) 



is induced on the sphere, and there is an accompanying 
image dipole in the sample 



Pirn = P 



(3.5) 



The electric field created by the image dipole at the cen- 
ter of the sphere is 

£ 1 g„ _ 1 r3 

^--^T^^8(^^-^ii^- ^'-'^ 

This is just the way the matrix elements of N should 
work; given an incoming field at the sphere, N generates 
the fields reflected from the sphere and scattered back to 
it by the sample. 

If the incident E field points in the z direction the 
field from the image dipole, which now is proportional to 
^Fo^io is twice as strong as in the previous, parallel case. 
Consequently we have. 



eT + 2 es + 1 4(fi + i?)3' 



(3.7) 



One can also evaluate the K integrals describing the 
coupling due to magnetooptic effects. The results can be 
written 



7^ = j^EE _ T^EE _ r^EE _ j^EE 
— ^10,11 — -"-11,10 ^ -'^ 10,1-1 ~ -"-1-1,10 



V2 



ei 



i?3 



eT — 1 

eT + 2 (es + l)2 8{d + Ry 



(3.8) 



The other elements of K vanish in this case. 

If an s polarized wave with wave vector q impinges on 
the system the wave incident on the dipole can be written 



™' = -2-/J[i + /'^(qii)]^ 



(3.9) 



where (f) denotes the azimuthal angle of q, the wave vector 
of the incident light. Because N is dia gonal it is simple 
to solve for the full vector a using Eq. ( 2.30 ), 



2T:i [Y r ... 



(3.10) 



The sphere is only polarized parallel to the surface at this 
stage {a[Q^ = 0). 

Before proceeding furthe r, let us insert the above result 
for a into Eqs. (2.3S) and (2.39) and calculate the polar- 



ization of the light sent out. Since E^J^ and E^J^ have 
many factors in common the final result can be written 
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PCP 



2Im 



S*P 



\S\' 



(3.11) 



where 



5" = 
P 



PpWw) 



and 



(3.12) 



With appropriate parameter values Eq. (3.11) predicts 
a rather large circular polarization even if the sample is 
non-magnetic. If pcp given by this expression is plotted 
as a function of (j)' one typically obtains an oscillating 
curve. We also note that in this limit, 4>' and (p' + tt 
yields exactly the same pcp- Thus, if we introduce the 
quantity 



Apcp(0') = Pcp(0') - Pcp{(p' + tt). 



(3.13) 



we find that in the dipole limit, without a Kerr effect, 
^PCP = 0. This is an important observation since we 
will see that when the sample is magnetic this symmetry 
is lost. 

Due to ei, there are two contributions to a'™*; the or- 
dinary Kerr effect upon reflection of the incident field 
yields a"^"^^ , 



-ei 



q\\P 



27r-\/3/47r sin </) 



Ps{p + Ps){esP + Ps) 



COSf 



and the tip Kerr effect yields 



Ka 



At:K sin 




Ps{q\\)] 



(3.14) 
(3.15) 



Both contributions are first order in ei, both have a sin;/) 
dependence on the direction of propagation of the inci- 
dent light, and the Kerr effect induces a dipole moment 
perpendicular to the sample surface. We also see that 
resonances in the tip-sample interaction described by the 
factor 1/(1 — iV||) can enhance the tip Kerr effect (though 
for a point dipole model for the tip these effects are usu- 
ally small). 

We can now solve for a' by making use of 



a' ={l-N 



Ka 



(3.16) 



where further tip-sample interactions are included. Next, 



inserting the sum a = a + a' into Eqs. (2.38) and (2.39) 



the polarization is found. 

Before analyzing this it is illustrative to see what re- 
sult we get for the magnetic contribution to the circular 
polarization, to lowest order in the ratio of E^^\ to E^^^^, 
and comparing this to our earlier result in the reverse 
situation with light being emitted by tunneling electrons 



instead.B Inserting the dipolar results above for the ordi- 
nary and the tip Kerr effect into Eqs. ( 2.3S ) and ( 2.39| ) 
we obtain to lowest order, retaining only magnetic con- 
tributions to the circular polarization 



p^i « Im 



x2 



sm sm 



1 + Gj 

1 + G|| \l + psiqii) 

2 Q/ 



cos (0 — (j)') ' 



-Pl(^) DHG\\ \ 1 + Pp(g||) 

(3.17) 



sm ( 



where Pl{^) is defined in Eq.(2.12) andlil 

2ei(co) 



(3.18) 



Furthermore the image factors G are defined through 

G||,^ = iV||,^/(l-7V|l,^) (3.19) 

for the parallel and perpendicular cases, respectively. 
The result for the circular polarization of light in the 
corresponding spontaneous emission STM configuration 



„STM 

Pcp 



2Im 



-Pl{^) , D{uj)G\\\ 1 



sin 6*' 



Pp 



(3.20) 



The structure is very similar to that of (3.17) but there 
are a couple of important differences. For the light-in 
situation the whole expression is basically larger by the 
factor G_L . This is one reason for why p^p can be several 
orders of magnitude larger than p^™ ; the perpendicular 
field enhancement, G_l, can be very large in the junc- 
tion between tip and sample, while G\\ is of order unity. 
Moreover, in the last factor (1 -t- ps) and (1 -f pp) have 
changed places. This is also an important reason for the 
large values reached by p^p- ^ metal surface (1 -I- p^) 
is usually considerably smaller than 1. 

Continuing the analysis of the full results we find that 
they can be summarized as 



{ia^ -f ay)/\/2 



(3.21) 



for a, where ax, ay and are proportional to the induced 
dipole moments in the x, y, and z directions, respectively. 
Let us assume from now on that the incident wave propa- 
gates in the negative y direction, so that cj) — —tt/2. This 
means that ax will be large, whereas ay and az are much 
smaller since they originate from the Kerr effect. There- 
fore looking at the radiation sent out in two opposite di- 
rections {(f)' and (j)' + tt) we see that the contributions to 

-^rad ^'^'^ -^rad causcd by the X and y components of the 
dipole moment both change sign upon changing the angle 
of observation, leaving the polarization state unchanged 
as we have already discussed. However, the new feature. 
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the dipole moment in the z direction adds the same con- 
tribution to E\J^ in both directions, thereby leading to 
an asymmetry in the angular dependence of the polar- 
ization. Anomalous reflection of the radiation sent out 
as a result of the dipole moment oscillating in the x di- 
rection also yields contributions to E\^^ and E^l^^ that 
do not change sign when changing 0' to ((>' + n. From 
this follows that the magnetism-related contribution to 
the circular polarization has a characteristic angular de- 
pendence which, as we will see, persists also when going 
beyond the dipole model. 



IV. RESULTS AND DISCUSSION 



A. Numerical results 



In this section we will present results of our numeri- 
cal calculations. We have solved the multiple-scatterin, 
problem at hand using the theory developed in Sec 
In these calculations we retained multipoles for which 
I < 30 and |m| < 3. This means that we get results that 
are essentially numerically exact; the relative accuracy 
obtained for Apcp is better than 10"^ when i? = 20 nm 
and d = 0.5 nm. If Imax is reduced to 10 the relative 
errors are about 1 %, and with Imax = 5 they amount to 
5-10 %. 

Figure ^ shows results for the polarization pcp as a 
function of the observation angle 0'. These results were 
obtained from a calculation using a silver tip and an Fe 
sample. For comparison, we have also performed a calcu- 
lation using a W tip [panel (c)] . The photon energy was 
1.6 eV and both the incident and the scattered light prop- 
agates in directions forming an angle of 1 radian with the 
surface normal. This means that 9 =1, while 9' = tt-I. 

Let us begin the discussion by looking at the results 
in Fig. 1^ (a), where the Ag tip radius was set to 10 nm. 
Results are shown both with (thick curve) and without 
(thin curve) the magnetooptic effects; the two curves are 
very similar indicating that "geometric effects" cause the 
dominating contribution to the polarization pcp- How- 
ever, a closer comparison reveals that the results includ- 
ing magnetooptic effects show a less symmetric angular 
dependence than those obtained without taking magne- 
tooptic effects into account. Therefore magnetooptic ef- 
fects can be isol ated more effectively if the quantity Apcp 
defined in Eq. ( 3.13| ) is plotted. We see that Apcp has 
peak values exceeding 10 %, and at the same time the 
polarization shows a simple, characteristic variation with 
the observation angle 0'. It appears that these effects are 
large enough to be detectable. We have also plotted the 
quantity Sp, i.e. the change in pcp due to magnetooptic 
effects. The two curves are very similar in shape. 

In Fig. ^ (b) the tip radius was increased to 20 nm. The 
tip is now large enough that phase differences between 
light scattered off different parts of the tip influence the 



polarization. This gives Apcp a more complicated angu- 
lar behavior, and it would be somewhat more difficult to 
pick out the part of the signal originating from the Kerr 
effect. For an even larger tip radius this effect becomes 
more pronounced. 

There are two different magnetism-related processes 
that contribute to the circular polarization of the scat- 
tered light. We already mentioned this in connection 
with Eqs. ( ^.32 ) and (3.17), and the same distinction 
was also made in our previous work on light emission. El 
To begin with there is an ordinary Kerr effect in which 
the polarization conversion, governed by the coefficient 
p'pg, takes place primarily when the incident wave first 
hits the sample surface. This induced p polarized wave 
is subsequently enhanced in the cavity between tip and 
sample; both tip and sample screen the electric field, sur- 
face charges are induced on both of them and this gen- 
erates secondary fields, etc. In the other process, which 
we have called the tip Kerr effect, the incident wave is 
reflected back and forth between tip and sample a num- 
ber of times and an electric fleld along the z axis is built 
up because one of the scattering events off the sample is 
anomalous. In this case the anomalously reflected wave 
typically is evanescent since it describes the near fleld 
around the tip. The conversion factor p' dominates 



these proces ses, b ecause, as can be seen from Eqs. ( 2.12 ), 



( p'.15 ), and ( 2.16 ), unhke Pp^ and p'^p it has a non-zero 
limit when 3> k. 

In the full (multipole) calculations the relative impor- 
tance of the two processes can be estimated, for example, 
by temporarily eliminating one or the other contribution. 
We find that the ordinary Kerr effect in the present case 
is responsible for some two thirds of the magnetic contri- 
butions to the circular polarization of the scattered light. 
In the calculation addressing emission of circularly polar- 
ized lighda we found the opposite situation; the tip Kerr 
effect dominated then. It is relatively straightforward 
to understand this difference qualitatively: Near-field ef- 
fects are considerably more important in the case of light 
emission because in that experiment one measures how 
much radiation is produced by a source, the tunnel cur- 
rent, that is basically localized to a nanometersized re- 
gion in space. In the present calculation we instead probe 
the degree to which light is scattered off the tip-sample 
system; the near-field effects play a role also here but it 
is not as pronounced. 

We have concentrated on studying the polarization 
properties of the scattered light at relatively low photon 
energies (< 2 eV) because ei for the magnetic materials 
takes the largest values there. In that case the results ob- 
tained with a W tip (shown in panel c) and with a Ag tip 
(panel a) are not too different from each other. This re- 
sult may at first seem surprising given that Ag is a much 
better conductor than W and has well-defined surface 
plasmon excitations. However, for the photon energies 
dealt with here the field enhancement does not occur as 
a result of truly resonant interactions with interface plas- 
mons. Instead the electric field in the tip-sample cavity is 
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enhanced by a less dramatic mechanism. Both Ag and W 
screen the electric field and therefore surface charges that 
interact with each other are built up on the tip and sam- 
ple. At higher photon energies a Ag tip certainly gives a 
stronger field enhancement than a W tip, but then again, 
since ei for the magnetic materials is much lower there, 
this frequency range is of limited interest in the present 
context. 

Figure ^ illustrates in more detail how the polariza- 
tion properties of the scattered light varies with photon 
energy. The highest degree of polarization is reached at 
1.5 eV. For photon energies beyond 2 eV the Kerr ef- 
fect has a rather small influence on the polarization, the 
main reason for this is that the magnitude of ei for Fe 
decreases with increasing hu. 

In Fig. ^we plot the differential scattering cross section 
for a Ag tip with a radius of either 10 nm or 20 nm. As 
a reference we have also plotted the result obtained for a 
Ag sphere with radius 10 nm in free space. The sphere 
in front of an Fe sample has a much smaller cross section 
than the one in free space because of destructive interfer- 
ence between the waves sent out directly from the sphere 
and those reflected from the sample, which in this case 
can be thought of as originating primarily from an image 
sphere. The radiation patterns differ also between the 
two cases. An isolated sphere displays a dipole pattern, 
but with the sample present, the angular distribution of 
radiation mainly follows a quadrupole pattern. To sec 
what the cross sections mean in terms of photon counts, 
let us assume that a laser with an intensity of 10^ W/cm^ 
is illuminating the tip. This is equivalent to a photon flow 
of ~ 10^ photons/(sec. nm^). In the case of the 10-nm- 
sphere in front of the Fe sample, a detector covering a 
solid angle of 0.1 steradian would collect ^ 10'^ photons 
per second. This intensity should be enough to perform 
a polarization analysis. 



B. Experimental implications 

The measurement of the magnetic contribution to the 
circular polarization should be feasible in the sense that 
the scattered light intensity is sufficient and the degree of 
circular polarization of the scattered light is large enough 
to be detectable with existing optical techniques. Fur- 
thermore, it should at least be possible to reach a resolu- 
tion of 10 nm. However, the removal of the background 
circular polarization which is not magnetic in origin is a 
nontrivial experimental problem. The Kerr contribution 
varies with scattering angle in a way which is different 
from other nonmagnetic contribution to the polarization, 
provided the tip is not too large. This should help in sep- 
arating the two contributions to the circular polarization. 
It may also be possible to use the fact that the Kerr con- 
tribution to the circular polarization is more sensitive to 
the tip-sample distance than is the background circular 
polarization, for example by letting the tip vibrate back 



and forth. To be specific, increasing the tip-sample dis- 
tance from 0.5 nm to 5 nm in Fig. ^ (a) decreases the Kerr 
contribution (Apcp) by about a factor of 2, whereas the 
background contribution is essentially unchanged. Vary- 
ing the photon energy could be yet another way of differ- 
entiating between Kerr and background contributions. 

In view of the large background, the size of the tip 
would be a major concern in designing an experiment in 
practice. If a traditional STM tip made from a metallic 
wire is used, regardless of how sharp the tip is, scattering 
will take place everywhere in the tip shaft. This scatter- 
ing will typically increase the geometric-background con- 
tribution to the circular polarization, and make it more 
difficult to detect the magnetic properties of the surface. 
Of course, the incident radiation would be focused to 
a relatively small part of the tip, but this may not be 
enough to obtain a high contrast between the contribu- 
tions to the polarization induced by the Kerr effect and 
geometric effects, respectively. 

As we have indicated in the Introduction, one possi- 
ble way of bypassing this difficulty could be to follow a 
scheme .closely related to that introduced by Silva and 
Schultzja using a small, isolated metallic particle on the 
surface of the magnetic sample as the probe. In order to 
scan different parts of the surface, the metallic particle 
would have to be moved around using an AFM. This kind 
of nanoparticle niaiiipulation has already been achieved 
in other contexts. E£l ■— ■ 

Pufall, Berger, and Schultzlij have measured the Kerr 
rotation of light scattered from Ag particles placed on a 
magnetic substrate. They found a Kerr rotation that had 
qualitative features (in terms of observation-angle depen- 
dence) in common with the results presented here. The 
measured Kerr rotation, however, was less than 1 mrad 
(corresponding to a magnetism-related degree of polar- 
ization of 0.1 % or less). It appears from Ref. |l9| that 
particles of a radius of at least 50 nm were used in the 
experiment and at the same time the incident laser light 
was tuned to the scattering resonance of the particles at 
a wavelength of 460 nm (w 2.7 eV). While both of these 
choices ensure that the scattering cross section is much 
larger than with the parameter values used in the calcu- 
lations above, this is achieved at the cost of a very small 
degree of magnetism-related polarization of the scattered 
light. First of all ei(w) in common magnetic metals such 
as Fe and Co decreases with increasing photon energy in 
the range of visible light (cf. Fig. ||). Secondly, with a 
larger particle the in-plane dipole becomes more effective 
in sending out radiation because the image sphere is fur- 
ther away from the real sphere. As a result the scattering 
cross section but also the nonmagnetic contributions to 
pcp increase dramatically. Indeed, using a silver parti- 
cle of radius i? = 50 nm on a Fe substrate and a pho- 
ton energy of 3 eV in our calculations we find that the 
magnetism-related contribution to the light polarization 
(i.e. Sp) is only « 0.1 %. With R =50 nm and a photon 
energy of 1.6 eV, (5p « 2 % at most, to be compared with 
6-7 % for i? = 10 nm (see Fig. | (a)). Thus, based on 
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the results found here we recommend that experiments 
of this kind should use smaller photon energies (1.5 eV 
or so) and smaller particles (~ 20 nm) than before. 
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APPENDIX: 

In this Appendix, we derive expressions for the func- 
tions / and g that were introduced in Sec. p^I D| . With a 
plane wave impinging on the sphere, the full electric field 
outside the sphere can be expanded as in Eq. ( p^ ), and 
the regular (ji) terms describe the incoming wave. We 
will calcu late the proportionality factors / appearing in 
Eq. (2.23), and therefore we first need to evalu ate t he a 
coefficients in Eq. ( ^.1^ ). It is clear from Eq. ( 2. IS ), in 
view of the orthogonality relations 



/ 



dfi X;* „j, (fi) • X;m(r2) = 6w6„ 



(Al) 



and 



dQX*,^,in) ■ V X [{ajiikr) + bhi{kr))Ximm] = 0, 

(A2) 

for the vector spherical harmonics, that 

kal^''jiikr) = Jdn (X„„)* • E. (A3) 

To calculate aj^ , it is easier to use the overlap between 
the B field and a vector spherical harmonic. 



(A4) 



To evaluate these surface integrals, we use the expansion 
of a scalar plane wave (i.e. the factor e''^ '' ) in terms of 
spherical harmonics and insert this in Eq. (2.1) yielding 



E(r) = {£;("'[z X q||] -£;(P)[qx (z x q||)]} e'l'^^p" 

oo 

x47r J2''ji{kR)Y,{-^rYi,-mmYi,m{^<i)- (A5) 



The incident B field can be expanded in a similar way. 
The factor e^'i'^^'p'^ compensates for the fact that we here 
use a coordinate system with the origin at the center of 



the sphere [Rsph = — z (R + d)]. The composite variables 
il and ilg denote the directions of f (i.e . the angles 9 and 
(f)) and q {9q and 0^). Equation (A5) is also valid for 
evanescent waves that decay exponentially in the posi- 
tive or negative z direction in which c ase co s Oq is purely 
imaginary. From the definition in Eq. ( 2.20 ) and the fact 
that L can be written in terms of ladder operators as 

L = i(x - iy)L+ + ^(x + iy)i- + zL^, 

it is clear that X/m can be expressed in terms of the 
spherical harmonics Y/.m+i, ^^/.m, and Yi „i-i- Combin- 
ing this with the expansion in Eq. (A5) yields the overlap 
integrals 

/,™"(q) = //™^(q) = k-' e'•'•''-^ (A6) 

and 

//™ (q) = - ftl'i^) = ViM e^'*-^-^ (A7) 
Here we have introduced 
27ri'(-l) 



C//m(q) 



N/^a+T) 



[^+F+{l,m)Yi,..m-iinq) + 



where the polarization vector 

^ = z X q|| y cos(?!)q 
C± = ± i^v = ±ie 
In a similar way 



with 



±i4 



(A8) 



(A9) 



27ri'(-l)' 



[T]+F+{l,m)Yi^^„i-i{^q) + 



via+T) 

+rj^F^{l,m)Yi^^^+i{nq) ~ rj,2mYi^^^{flq)] , (AlO) 

where the polarization vector 

fj ~ fc^-^q X (z X q||) 

— x(— cos 9q COS (j)q) + y(— COS 9q sin i/ig) -I- z sin 9q, with 
r]±=r]^±iTjy ^ - cos 9qe^''''''. (All) 



In both Eq. (|A8[ ) and Eq. (|A10D , F+ and F_ are short- 
hand symbols for the numerical factors produced by the 
ladder operators and L_, thus 



F+{1, m) — {I — m){l -|- m -I- 1), and 
F_{l,m) = y/{l + m){l-m + l). (A12) 

Finally, we note that the (jjq dependence of both and 

Vim follows 



C/h„(q) oc e-""^^ V,,„(q)(xe-™^''. 



(A13) 



The same is true for fi„i (provided Rsph lies on the z 
axis) . 
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Next we must deal with the overlap going in the other 
direction, from a given multipole to a plane wave. The 
vector equivalents of the Kirchoff integrals provides one 
way of doing this. They read 

E(r) = E'=''*(r) + J dS' [ikc{n' x B(r'))G(r, r') 

+ (n' X E(r')) X V'G + (n' • E(r'))V'G] (A14) 

and 



B(r) =B°'^'(r) 



dS' 



--(n' X E(r'))G(r,r') 



+ (n' X B(r')) X V'G + (n' • B(r'))V'G] . (A15) 

Here G denotes the Green's function of the scalar 
Helmholtz equation in free space, thus G solves 



[V2 + fc2]G(r,r') = -(5(3)(r_r'), 
where k = lo/c, and can be written 

G(r,r')- 



47r r — r' 



(A16) 



(A17) 



To get t he tot al field at a point in free space, the integrals 
in Eqs. ( A14 ) and ( |A15| ) should be over the surfaces of 
all the scatterers that are present, and the fields entering 
the integrals should be the exact fields at those interfaces. 
Here we restrict the attention to the fields scattered from 
the sphe re. I n that case , only the Hankel function terms 
in Eqs. ( 2.18 ) and ( 2.19 ) contribute. Evaluating the sur- 
face integrals using an expansion of the Green's function 
in terms of plane waves (propagating and evanescent). 



G(r,r') 



»qir(ri|-r;i) 



(A18) 



one finds that the field outside the sphere can be written 
d% 



lm,(T 



1 

(27r)2 



X [gf^(q)(z X q||) - C(q)(q X X qii))] e^'^"- (ai9) 

Here cr denotes a polarization [(E) or (M)]. The coupling 
factors are found to be 



5r™(q)=.9fr(q) 



-1 



d+m+l 



iq-Rsph 



(A20) 



and 



fff™(q) 



where Rgph is a vector pointing to the center of the 
sphere. As for the (pq dependence of gim, it is in view 
of Eq. (|A13D clear that 
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FIG. 1. The model geometry used in the calculation. The vectors q', p, and s indicate propagation and polarization 
directions in the radiation zone. 6 and (f) and 9' and 0' are the angles of incidence and observation, respectively, in a spherical 
coordinate system. 




M M ^^P m(S) ^pp 

FIG. 2. Illustration of the different "anomalous" scattering processes involving the Kerr effect that can take place at the 
sample surface. Conversion of s to p polarization and vice versa is most effective when the in-plane component of the wave 
vector is coUinear with the sample magnetization M. However, the Kerr contribution to the p to p reflection is most effective 
when qy and M are perpendicular (NB. In the figure to the right the viewpoint, but not the magnetization differs from the 
other figures). 



13 




p' p 



FIG. 3. Illustration of the general scheme employed in the calculation. An incident plane wave can be expressed as a 
multipole expansion of waves impinging on the sphere, either directly (o^"^) or after an initial reflection off the sample (a"''). 
Subsequently, multiple scattering takes place between the sphere and the sample. The sphere response functions s determines 
the strength of the outgoing spherical waves bi'm', which can be expanded into plane waves according to the functions gi'm' 
and after reflection off the sample (governed by p), the functions fim determine the strengths ai,n of the different multipole 
contributions to the reflected wave that again impinges on the sphere. Anomalous (Kerr effect) reflections which are governed 
by the functions p' generate waves incident on the sphere with strength a[^. 




14 



Polarization (%) 

Fe sample- Ag tip (R=20 nm) 
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FIG. 4. Calculated degree of polarization for a Ag tip probing a Fe sample. The curves in panel (a) were calculated with 
a tip radius -R = 10 nm, the ones in panel (b) with R — 20 nm. For comparison results calculated with a W tip with R — 
10 nm are displayed in panel (c). In all cases the tip-sample separation d —0.5 nm. The photon energy was set to 1.6 eV, 
(j) = —90°, and the polar angles of incidence and observation 8 and 6' , respectively, both lie 1 radian off the surface normal (i.e. 
= 1 rad. while n — 6' = 1 rad.). The two curves with the highest amplitude show the results for pcp- The thinner of these 
curves was calculated without accounting for the magnetooptic properties of the sample, whereas they have b een in cluded in 
the calculation yielding the thicker curve. The remaining curves display (i) the quantity Apcp defined in Eq. ( j.l3 ), and (ii), 
in panel (a), the difference 5p between the results with and without magnetooptic effects. 
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FIG. 5. The polarization asymmetry Apcp, defined in Eq. (3.13), plotted as a function of observation angle 0', for a rmmber 



of photon energies for a Fe sample scanned by a Ag tip with a radius of 10 nm, and with d set to 0.5 nm. 
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FIG. 6. Calculated differential scattering cross section for a Ag sphere in front of a Fe sample (full curves), and in free 
space (dotted curve) as a function of azimuthal observation angle. The photon energy is 1.6 eV, d — 0.5 nm, the sphere radius 
is 10 nm (lower full curve and dotted curve), and 20 nm (upper full curve), respectively. The polar angles of incidence and 
observation and 9' , respectively, both lie 1 radian off the surface normal (i.e. 9 = 1 rad. while n — 9' = 1 rad.). 
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TABLE I. The ofF-diagonal element ei(t<j) of the dielectric 
tensor for Fe and Co for a few photon energies. The data are 
taken from Ref. 



hiu (cV) 


Fe 


Co 


1.5 


-l.42-0.089i 


-1. 16-0.0891 


2.0 


-0.67-0.26i 


-0.40-0.151 


2.5 


-0.26-0.22i 


-0.20-0.0961 


3.0 


-0.13-0.12i 


-0.10-0.0671 



